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We know that the great circle is a circle whose plane passes through the centre of sphere. The part of a great circle 
on a sphere is known as a great circle arc. The length of minor great circle arc (i.e. less than half great circle) 
joining any two arbitrary points on a sphere of finite radius is the minimum distance between those points. Here 
we are interested in finding out the minimum distance or great 
circle distance between any two arbitrary points on a spherical 
surface of finite radius (like globe) for the given values of 
latitudes & longitudes. 


Let there be any two arbitrary points A(@, ,A,) & B(@2 , Az) 
lying on the surface of sphere of radius R & centre at the point 
O. The angles of latitude @, & @2 are measured from the 
equator plane & the angles of longitude A, &A, are 
measured from a reference plane OPQ in the anticlockwise 
direction. Here, we are to find out the length of great circle 
arc AB joining the given points A & B. Draw the great circle 
arcs AD and BC passing through the points A & B 
respectively which intersect each other at the peak (pole) 
point P & intersect the equatorial line orthogonally (at 90°) at 
the points D & C respectively. (As shown by the dashed arcs 
PD & PC in the figure 1) 


Join the points A, B, C & D by the dashed straight lines 
through the interior of sphere to get a plane quadrilateral 
ABCD and great circle arc BD. 


Figure 1: The dashed great circle arcs PD and PC passing through 
Now, the angle between the orthogonal great circle arcs BC two given points A & B, intersecting each other at the peak (pole) 


& CD, subtending the angles a=$, & B=A,—A, respectively on a spherical surface of finite radius R. 
respectively at the centre O of the sphere, meeting each other 

at the common end point C, is 8 = / . Now, consider the 

tetrahedron OBCD formed by joining the points B, C and D to the centre O (see the fig-1). A diahedral angle is 

the angle between two intersecting planes measured in a plane perpendicular to the both the intersecting planes. 

Now, the diahedral angle say 8 between the lateral triangular faces BOC and COD is given by HCR’s Inverse 

Cosine Formula [1] as follows 


cos a’ — cosacosB 
cos 8 = ——_—____— jintninatevade(L) 
sina sinB 
where, is 8 is the diahedral angle between lateral triangular faces BOC and COD intersecting each other at the 
line OC which is equal to the angle between great circle arcs BC and CD intersecting each other perpendicularly 


ie. 9 = = which is opposite to a’ = ZBOD and 
a’ =ZBOD, a= ZBOC, B = ZCOD are the vertex angles of triangular faces BOD, BOC and COD respectively 


meeting at the apex O of tetrahedron OBCD and angle a’ is always opposite to the diahedral angle 6 (as shown in 
the above fig-1 and fig-2 below). 
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point P, meet the equator orthogonally at the points D & C 


Now, substituting the corresponding values i.e. 0 = a= 
ZBOD,a = ZBOC = @¢2 and B = ZCOD = A, — A, in the 
above Eq(1) (HCR’s Inverse Cosine Formula) as follows 

m cosZBOD — cos #2 cos(Az — A,) 


oa. sin dz sin(Az — A,) 


cos ZBOD — cos ¢, cos(Az — A,) = 0 


cos ZBOD = cos 2 cos(Az — A,) ..(2) 


Similarly, from the figure-2, the diahedral angle say y 
between the lateral triangular faces BOD and COD of 
tetrahedron OBCD is obtained by _ substituting the 
corresponding values i.e. 9 = y which is opposite to a’ = 
ZBOC = $2, a = ZBOD and 6 = ZCOD = Az — A, in the 
above Eq(1) ( HCR’s Inverse Cosine Formula) as follows 


cos @, — cos ZBOD cos(A, — A,) 


Figure 2: The diahedral angle between the triangular faces BOC 
and COD of tetrahedron OBCD is equal to the angle between 
great arcs BC and CD which intersect each other at right angle 
because the great arc BC intersects equatorial line (i.e. great arc 
CD) at right angle. 


er sin ZBOD sin(A, — A,) 
cos @, — cos d@, cos(A, — A,) cos(A, — A,) ; 
= Sett ] f ZBOD fi Eq(2 
sin ZBOD sin, — 2,) (Setting value of cos rom Eq(2)) 
_ cos 2 (1 —cos*(Az—44)) _ cos@zsin*(Az — 1) _ cos 2 sin(Az — A;) 
~  sinZBODsin(A,—”,) — sinZBODsin(A, —”,) _ sin ZBOD 
cos @2 sin(a, — A 
cos y = cos b2 sin(Az ~ Ai) > (3) 


sin ZBOD 


Now, from the figure-3, consider the tetrahedron OABD formed 
by joining the points A, B and D to the centre O (also shown in 
the above fig-1). Now, the diahedral angle say 5 between the 
lateral triangular faces AOD and BOD is obtained by substituting 
the corresponding values i.e. 9 = 5 which is opposite to @’ 
ZAOB, a = ZAOD = ¢, and B = ZBOD in the above Eq(1) ( 


HCR’s Inverse Cosine Formula) as follows 


__ cos ZAOB — cos ¢; cos 2BOD 


sd sin @; sin ZBOD 
(90° = cos ZAOB — cos ¢, cos ZBOD 
a om sind, sin ZBOD 
(ey +6 = 90°) 


Sf 

ws 
Diahedral 

angle 


5=90°-y 


Figure 3: The diahedral angle 6 between the triangular faces 
AOD and BOD is equal to the angle between great arcs AD 
and BD which is opposite to the face angle 2AOB of 
tetrahedron OABD. 
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; cos ZAOB — cos ¢, cos 2 cos(A, — A,) (Setti eat aon Ba(2)) 
s--_eo(—____—e—e—e—e—————————————————————— ee eee Z 
siny sin @, sin ZBOD etting value of cos rom Eq(2) 


sin fd, siny sin ZBOD = cos ZAOB — cos ¢, cos #2 cos(Az — A,) 
cos ZAOB = cos ¢, cos @2 cos(A, — A,) + sing, siny sin ZBOD 
= cos ¢; cos @2 cos(A, — A,) + sing; sin ZBOD sin? y (-O0<y<zm) 


= cos ¢, cos d2 cos(A, — A,) + sing, sin ZBOD ./1 —cos?y 


= (A, —a,) + sing, sinzBoOD |1 £08 2 sin(Az — ax)" Cemne) 
= cos Pd; cos Pz cos(Az 1) + sind, sin in ZHOD rom Eq(3) 


Gaines upon sin? ZBOD — cos? #2 sin? (Az — a4) 
= — ZL. Ss 
cos P; cos Pz cos(Az 1) + sind, sin sin? ZBOD 


= cos d, cos Pz cos(Az — A) + sind, J1 — cos? ZBOD — cos? ¢2 sin? (Az — a4) 
Substituting the value of cos 2BOD from the above Eq(2), 


cos ZAOB = cos ¢, cos @, cos(A, — A,) + sing, 1 — (cos dz cos(A, — A,))? — cos? d, sin?(A, — A,) 


= cos ¢, cos dz cos(Az — A,) + sin dy ft — cos? d, cos*(Az — A,) — cos? 2 sin? (Az — a4) 


= cos g, cos @, cos(A, — A,) + sing, 1 — cos? fp (sin? (A, — A,) + cos?(A, — A,)) 


= cos ¢, cos @2 cos(A, — Az) + sing, V1 — cos? 2 (1) 
= cos d, cos P2 cos(A, — A,) + sin b, J sin? 2 
= cos g, cos @, cos(A, — A,) + sing, sin dz (“Os 


=> ZAOB = cos"1(cos ¢, cos g, cos(A, — A,) + sing, sing) .. (4) 


The above ZAOB is the angle subtended at the centre O by the great circle arc AB joining the given points A and 
B lying on the sphere (as shown in the above fig-1). Therefore the minimum distance between two points A and 


B on the sphere will be equal to the great circle arc AB given as follows 
The length of great circle arc AB = Central angle x Radius of sphere = ZAOB XR 
= cos /(cos ; cos fz cos(A, — A,) + sing, sing,) x R 


= Rcos (cos ¢; cos @2 cos(Az — A,) + sin d, sin dz) 
V O< 1,62, |A2 —Ay| <7 


The above formula is called Great-circle distance formula which gives the minimum distance between any two 


arbitrary points lying on the sphere for given latitudes and longitudes. 


Copyright© H.C. Rajpoot 


NOTE: It’s worth noticing that the above formula of Great-circle distance has symmetry i.e. if 6, & dz and 
A, & Az are interchanged, the formula remains unchanged. It also implies that if the locations of two points for 
given values of latitude & longitude are interchanged, the distance between them does not change at all. 


Since the equator plane divides the sphere into two equal hemispheres hence the above formula is applicable to 
find out the minimum distance between any two arbitrary points lying on any of two hemispheres. So for the 
convenience, the equator plane of the sphere should be taken in such a way that the given points lie on one of the 
two hemispheres resulting from division of sphere by the reference equator plane. 


Since the maximum value of cos~+(x) is 1 hence the max. of min. distance between two points on a sphere is 


= R(t) = 1R = half of the perimeter of a great circle passing through given points 


Case 1: If both the given points lie on the equator of the sphere then substituting ¢, = $2 = 0 in the great-circle 
distance formula, we obtain 


din = R cos *(cos 0 cos 0 cos(Az — A,) + sinOsin0) = Rcos *(cos(A, —A,)) = RA, —A,| = RIAA 


The above result shows that the minimum distance between any two points lying on the equator of the sphere 
depends only on the difference of longitudes of two given points & the radius of the sphere. In this case, the 
minimum distance between such two points is simply the product of radius R and the angle AA between them 
measured on the equatorial plane of sphere. This 


If both the given points lie diametrically opposite on the equator of the sphere then substituting |AA| = 7 in above 
expression, the minimum distance between such points 


R|AA| = R(t) = 7R = half of the perimeter of a great circle passing through thegiven points 


Case 2: If both the given points lie on a great circle arc normal to the equator of the sphere then substituting A, — 
A, = AA = 0 in the formula of great-circle distance, we obtain 


dinin = R cos~1(cos ¢, cos f2 cos(0) + sin d, sin dz) 
= Rcos"}(cos @; cos 2 + sind, sing.) = Rcos~1(cos(¢, — $2)) = Rid, — $2! 


In this case, the minimum distance between such two points is simply the product of radius R and the angle 
|@1 — @2| between them measured on plane normal to the equatorial plane of sphere. 


Case 3: If both the given points lie at the same latitude of the sphere then substituting ¢; = @2 = ¢ in the formula 
of great-circle distance, we obtain 


dinin = R cos! (cos f cos f cos(A, — A,) + sing sin @) = Rcos~1(cos? @ cos(A, — A,) + sin? d) 


In this case, the minimum distance between such two points is dependent on the radius and both the latitude, and 
longitudes. 
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Illustrative Example 


Consider any two arbitrary points A & B having respective angles of latitude @,; = 40° & 2 = 75° & the 
difference of angles of longitude AA = Az — Ay = 55° on a sphere of radius 25 cm. The minimum distance 
between given points lying on the sphere is obtained by substituting the corresponding values in the above great- 
circle distance formula as follows 


dmin = R cos~*(cos d, cos @2 cos(A, — A,) + sind, sin hz) 
= 25 cos +(cos 40° cos 75° cos(55°) + sin 40° sin 75°) 
= 25cos_1(0.7346063699582707) = 18.64274952833712 cm 


The above result also shows that the points A & B divide the perimeter = 27(25) ~ 157.07963267948966cm 
of the great circle in two great circles arcs (one is minor arc AB of length = 18.64274952833712 cm & other 
is major arc AB of length = 138.43688315115253 cm) into a ratio = 


18.64274992893712 1, 35 43688315115253 © 1:7-4 


Conclusion: It can be concluded that the analytic formula of great-circle distance derived here directly gives the 
correct values of the great-circle distance between any two arbitrary points on the sphere because there is no 
approximation in the formula. This is extremely useful formula to compute the minimum distance between any 
two arbitrary points lying on a sphere of finite radius which is equally applicable in global positioning system. 
This formula is the most general formula to calculate the geographical distance between any two points on the 
globe for the given latitudes & longitudes. This is a high precision formula which gives the correct values for all 
the distances on the tiny sphere as well as the large spheres such as Earth, and other giant planets assuming them 
the perfect spheres. 


Note: Above articles had been derived & illustrated by Mr H.C. Rajpoot (B Tech, Mechanical Engineering) 
M.M.M. University of Technology, Gorakhpur-273010 (UP) India Aug, 2016 
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